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ABSTRACT 

This paper presents a sensitive and comprehensive IRAC 3 — 8 yum photo- 
metric survey of white dwarfs for companions in the planetary mass regime with 
temperatures cooler than the known T dwarfs. The search focuses on descen- 
dents of intermediate mass stars with M > 3 M© whose inner, few hundred AU 
regions cannot be probed effectively for massive planets and brown dwarfs by any 
alternative existing method. Furthermore, examination for mid-infrared excess 
explores an extensive range of orbital semimajor axes, including the intermediate 
5 — 50 AU range poorly covered and incompletely accessible by other techniques 
at main sequence or evolved stars. Three samples of white dwarfs are chosen 
which together represent relatively young as well as older populations of stars: 
9 open cluster white dwarfs, 22 high mass field white dwarfs, and 17 metal- rich 
field white dwarfs. In particular, these targets include: 7 Hyads and 4 field white 
dwarfs of similar age; 1 Pleiad and 19 field white dwarfs of similar age; van 
Maanen 2 and 16 similarly metal-rich white dwarfs with ages between 1 and 7 
Gyr. 

No substellar companion candidates were identified at any star. By demand- 
ing a 15% minimum photometric excess at 4.5 /im to indicate a companion de- 
tection, upper limits in the planetary mass regime are established at 34 of the 
sample white dwarfs, 20 of which have limits below 10 Mj according to substel- 
lar cooling models. Specifically, limits below the minimum mass for deuterium 
burning are established at all Pleiades and Hyades white dwarfs, as well as simi- 
larly young field white dwarfs, half a dozen of which receive limits at or below 5 
Mj. Two IRAC epochs at vMa 2 rule out T > 200 K proper motion companions 
within 1200 AU. 
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1. INTRODUCTION 

The first strong candidate and certain substellar objects identified outside the Solar 
system were all discovered orbiting e volved degenerate stars: the probable brown dwarf cora - 
panion to the white dwarf GD 165 (IKirkpatrick et al.lll999l : iBecklin. fc Zuckerman.1 Il988l) . 
and the planetary system orbiting the pulsar PSR 1257+12 (jWolszczan fc Fraillll992l ). Fur- 
thermore, the first directly detected, unambiguous substellar and planetary mass objects 
were imaged as wide companio ns orbiting low lumin osity primaries: the brown dwarf sec- 
ondary to the M dwarf Gl 229 (INakajima et al.l 119951), and the planetary mass secondary t o 
the young brown dwarf commonly known as 2M1207 (ISong et al.ll2006l : IChauvin et al.ll2005l ). 
Historically, as well as astrophysically, these properties provide a clear advantage over other 
types of primaries in the quest to directly detect radiation from bound substellar objects of 
the lowest mass, such as planets. 

Observations indicate that bound substellar objects and planetary system components 



cent giant stars known to harbor substellar and planetary companions ( 


Nied 


zielski et al. 


2007; 


Reffert et al. 


2006: 


Hatzes et al. 


2006. 


2005; 


Sato et al. 


2003; 


Frink et al. 


2002 


). Sec- 



ond, there a re at least 3 white dwarfs which have close, unaltered and unevolved substella r 
companions (IBurleigh et al.ll2006l : iMaxted et al.ll2006l : iFarihi. Becklin. fc Zuckermarul2005al ). 
Third, there are at least 10 white dwarfs with infrared excess due to debris disks which in- 



dicate a growing probability of orbiting rocky planetesimals (IFarihi. Zuckerman. fc Becklin 



200a lJura et al.ll2007bl : I Jura. Farihi. fc Zuckermanll2007al : iKihc fc Redfieldll2007l : iReach et al. 



20051 ). These cool white dwarfs with warm orbiting dust also display anomalous photo- 
spheric metals which are almost certainly accreted from their circumstellar material. The 
origin of the orbiting material and the dynamical interactions necessary to bring it close 
enough to the star to be accre ted, are consistent with remnant planetary systems (jjura 
2OO3I : bebes fc Sigurdsso^l2002h . 



The earth-size radii and consequent low luminosities typical of white dwarfs are clear 
advantages when searching for light emitted from cold jupiter-size planets and brown dwarfs. 
Prior to the launch of Spitzer, generally speaking, only M and L dwarf companions could be 
detected directly as excess infrared emission from white dwarfs (ITrembley fc Bergeron! 120071 ; 
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Hoarc 


et al. 


2007; 


Burleieh et aliboOC 


i: Farihi. Hoard. & Wachtei 


2006; 


Farihi. Becklin. & Zuckerman 


2005a 


Farihi. Zuckerman. & Becklin 


2005b; 


Farihl2004: 


Wachter et al. 


2003 


Zuckerman & Becklin 




1987 




Probst 


1983). As a benchmark, a typical 10000 K degenerate and an L5 dwarf 





are ab out equally luminous at fC-band (jPahn et al.l |2002| ; iBergeron. Saumon. fc Wesemael 



1995bj), while ground-based observations of white dwarfs at longer wavelengths (where the 



contr ast for coole r companions would improve) are prohibited by overwhelming sky bright- 



ness (lGlasslll999l ). 



Owing to the capabilities of Spitzer (IWerner et al.ll2004l ). a Cycle 1 IRAC ( IFazio et al. 



20041 ) program was undertaken to photometrically search for massive planets and cold sub- 
stellar companions to relatively young and old white dwarfs, respectively. Specifically, the 
target sample includes white dwarfs in the Hyades and Pleiades, high mass field white dwarfs, 
and metal-rich field white dwarfs. Farihi, Zuckerman, & Becklin (2008; Paper I) describe 
photometry for all the older (metal-rich) degenerate targets, while this paper presents a 
synopsis of the IRAC results for the younger (open cluster and high mass) degenerates, and 
upper mass limits for unresolved companions to all the Cycle 1 targets. 



2. RELATIVELY YOUNG AND OLD DEGENERATE TARGETS 



While highly evolved, white dwarfs are not necessarily old. This fact is exemplified by 
the nearest and brightest degen erate star in the sky, Sirius B, with a mass of M ^ 1.00 Mq 
and a total age of t ^ 240 Myr (ILiebert et al.ll2005bl ). For stars which evolved essentially as 



mass, derived primarily from studies of open clusters ( 


Dobbie et al. 


2006a; 


Kalirai et al. 


2005; 


Williams, Bolte. & Koester 


2004; 


Claver et al. 


2001; 


Weidemanr 


2000 


, |l990. 


1987 


)• 



This initial-to- final mass relation yields an estimate of the total age for any particular white 
dwarf provided its mass and cooling age are accurately known. For white dwarf masses below 
0.6 Mq, the initial-to-final mass relation is quite s teep and small errors in degenerate mass 



can lead to large errors in main sequence lifetime (IBurleigh. Clarke, fc Hodgkinll2002l ). 



2.1. Hyad and Pleiad Targets 



The Hyades and Pleiad es are relatively 



(i = 46 pc and r = 625 Myr (IPerrvman et al 



young 



1998 



and nearby open cluster s; the former at 



at r = 125 Myr and d = 13 2 pc ( Stauffer. Schultz. &: Kirkpatrickl 



Pinsonneault et al.l 1998), and the latter 



1998 



Pinsonneault et al. 



19981 ; ISoderblom et al.lll998l ). Classically, the only white dwarf Pleiad is EG 25 (LB 1497), 
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but recently the possibility has been raised that the massive white dwarfs G D 50 and PG 
0136+ 251 may have originated in the same region which gave rise to the cluster (jPobbie et al. 



2006bl ). The latter stars are considered in the following section, and only EG 25 is listed 



among cluster targets. Table [T] lists all observed open cluster white dwarfs, including the 
seven classical single white dwarf Hyads and EG 265 (y411 r), which is either a proper 
cluster member or part of the Hyades supercluster (lReidlll993l : lEggenlll984l ). Excluded are 
the white dwarf Hyads in binary systems; V471 r and HZ 9. The open cluster targets come 
from this study {Spitzer Program 3549; PI Becklin), with the exception of Hyades targets EG 
39 and EG 42, which were extracted from the Spitzer archive (Program 2313; PI Kuchner). 



2.2. High Mass Targets 



As with Sirius B and the white dwarf Pleiad, a white dwarf mass near 1.0 Mq implies 
a short main sequence hfetime for unadulterated single star evolution, regardless of asso- 
ciation. More specifically, were these young white dwarf s identical in mass (depending on 



the r eference, their masses differ by no more than 10%; iLiebert et al.ll2005bl : IClaver et al. 



200ll ). their total age difference can be estimated from the difference in their effective tem- 
peratures; Ted = 25, 200 K for Sirius B, and Tes = 31, 700 K for EG 25. For these temper- 
atures, log g = 8.6 (very nearly 1.0 Mq) hydrogen atmosphere models predict cooling ages 



difference ( 


Barstow et al. 


2005; 


Claver et al.l 


2001; 


Bergeron, Saumon, & Wesemael 


1995b: 


Bergeron, Wesemael, & Beauchamp 


1995c: 


Bergeron, Liebert, & FuUbright 


1995a 


). In real- 



ity the total difference in their ages is due both to differential cooling and unequal main 
sequence progenitor lifetimes. 

Any hot, high mass white dwarf which evolved as a single star will be similarly young, 
or even younger for higher temperatures or masses. In this paper, it is assumed that all 
high mass white dwarfs are descended from single main sequence star progenitors of in- 
termediate mass, but this may not be the case. There appears to be indirect evidence in 
favor of, as well as against, th e existence of high mass white dwarfs resulting from mergers , 



but no direct evidence exists (Hansen, Kulkarni, fc Wiktorowiczl l2006l: iFerrario et al.l 12005 



Liebert, Bergeron, fc Holberg]|2005al ). Ten hot and massive white dwarf targets come from 
this study {Spitzer Program 3549; PI Becklin) and one dozen similar targets were extracted 
from the Spitzer archive (Program 3309; PI Hansen). Table [2] lists the 22 hot field white 
dwarfs with masses M > 0.9 Mq selected for study. 
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2.3. Older Targets 

Included in the sample are metal-rich white dwarfs from Paper I [Spitzer PTogi am 3548; 
PI Zuckerman) with the addition of vMa 2, which was extracted from the Spitzer archive 
(Program 33; PI Fazio). These targets are relatively old, with total ages between 1 — 7 
Gyr and are listed in Table [3] for completeness, although their IRAC fluxes are previously 
published in Paper I with the exception of vMa 2. 



3. OBSERVATIONS AND DATA 



For white dwarf targets in Programs 3548 and 3549, the details of the IRAC observing 
strategy, data reduction and analysis are described in full detail in Paper I. In these programs 
a total exposure time of 600 s was utilized for each target in all bandpasses. For white dwarf 
targets extracted from t he Spitzer archive, the exposure times were shorter and occasi onally 
in only 2 bandpasses (see iMuUally et al.ll2007l : lHansen. Kulkarni. fc Wiktorowiczll2006l ). For- 
tunately, all targets were unambiguously detected at 4.5 /im, the wavelength which places the 
best constraints on spatially unresolved cold substellar and ma ssive planetary companions, 
according to models for th e appropriate range of ages and masses ( jBurrows. Sudarsky. fc Lunine 
2OO3I : ISaraffe et al-lboosh . 



3.1. Photometry and Upper Limits 

Paper I contains a detailed discussion regarding the consistency of measured IRAC fluxes 
of white dwarfs compared to model predictions and concludes the photometric accuracy is 
well described as 5%. Also described there is a conservative approach to estimating the 
signal-to-noise of IRAC detections in the presence of possible confusion and spatially varying 
background. In Paper I, all targets were detected at all wavelengths, which is not the case 
here. 

To create upper limits for nondetections, aperture photometry was performed as de- 
scribed in Paper I at the nominal location of the white dwarf, derived from one or more 
IRAC channels in which the source was positively detected. Utilizing the smallest radius 
(r = 2 pixels) for which there are published aperture corrections, the flux in this aperture 
was compared to the per pixel sky noise multiplied by the area of the aperture, and the larger 
of these values was taken to be the upper limit, after an appropriate aperture correction. 
In nearly all cases, the larger value was given by the additive noise in the aperture, but 
there were a few cases in which there was flux measured above this level. In these cases. 
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it appeared possible or likely that the measured flux originated from background sources 
as evidenced by pixel shifts in the centroid of the source flux compared to the other IRAC 
channels, or sources which were apparently extended. 

It is noteworthy that the photometric errors used here, as well as the derived upper limits 
for no ndet ections, are soniewhat larger than those published in lHansen. Kulkarni. fc Wiktorowic2 
(120061 ) and lMuUally et al.l (120071 ) for the same observations. In some cases where white dwarf 
flux is reported by those authors. Table S] indicates only an upper limit for the reasons stated 
at the end of the previous paragraph. IRAC fluxes and upper limits for 32 of 48 studied 
stars (excepting those previously published in Paper I) are listed in Table H] and plotted in 
Figures [THU 



4. ANALYSIS 
4.1. Total Ages 

To infer companion mass limits using substellar cooling models, each white dwarf target 
requires an assessment of its total age since formation as a main sequence star. Substellar 
companions which form in a binary will be truly coeval, while massive planets are thought 
to form within 10 Myr of their main sequence hosts. For the open cluster white dwarfs, their 
total age is the cluster age, while for field white dwarfs the following methods are used to 
estimate their ages. 

First, any field white dwarf with a mass and effective temperature similar to or higher 
than the Pleiad EG 25 is assumed to be of similar age, roughly 0.1 Gyr. There are sufficient 
uncertainties in both the masses and temperatures of these stars, as evidenced by the 10—20% 
variation among parameters cited in the literature for the same objects (see Table [2] for a 
list of references) which translates into errors in cooling ages on the order of 10 — 20 Myr. 
Fortunately, this type of error should be offset when assessing a total age because higher mass 
white dwarfs cool more slowly (smaller surface area) yet have shorter inferred main sequence 
lifetimes than their less massive counterparts. For the high mass field stars considered here, 
the modest errors in cooling age and inferred main sequence lifetime, which result from 
uncertainties in white dwarf parameters, are comparable in magnitude and therefore tend 
to cancel out. Given these uncertainties for the hot and massive field white dwarfs, it seems 
prudent to assign a 20% uncertainty in their total ages, or r = 0.125 ± 0.025 Gyr. 



Second, all cooler and less massive (i.e. older) fie ld degenerates have their tot a l ages 
assessed following the procedure employed originally bvlBurleigh. Clarke, fc HodgkinI (120021 ) 
and more recently by lDebes. Sigurdsson. fc Hansen! (120071 ). This latter method utilizes the 
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initial-to- final mass to relation to obtain a main sequence mass from the current, known 
white dwarf mass. A main sequence lifetime is then assigned based on the inferred main 
sequence mass, and this is added to the white dwarf cooling lifetime to obtain an approximate 
total age 



wd 



Cooling ages come from models of P. Bergeron (200 2, private communication; iBergeron. Saumon. fc Wesem; 

1995bl : iBergeron. Wesemael. &: Beauchamdll995cf). while m ain sequence lifetimes were cal- 

culcated using the formulae of iHurley. Pols. &: ToutI (120001 ). Tables [5] and [6] list the relevant 

ages for all targets together with upper mass limits for substellar companions determined 

as described below. It should be stated that this general method is the best available to 

estimate the total age of white dwarfs not belonging to open clusters or multiple systems 

from which another age constraint might be gleaned. However, there are several sources 

of uncertainty in the estimation of total ages, including but not limited to: the slope of 

the initial-to-final mass relation, assumed main sequence lifetimes, and white dwarf model 



uncertainties; specifically, spectroscopic parameter fits and cooling ages (jKalirai et al.l 12008 



Dobbie et al.ll2006al : iFerrario et al.ll2005l ). Owing to these facts, the uncertainty in the total 
ages of older field white dwarfs is taken to be 25%. 



4.2. Unresolved Companion Mass Limits from 4.5 fim Fluxes 

All white dwarfs in Tables [THS] have measured IRA C fiuxes at 4.5 /zm, where cold (Tpff < 



1000 K) substellar objects are predicted to be brightest (IBaraffe et al.ll2003l : lBurrows. Sudarsky. fc Lunine 



20031 ). This bandpass is best for placing limits on any spatially unresolved substellar com- 
panions. The fiux errors in Table H] are la values, but for reasons discussed in Paper I, and 
to be conservative, an unambiguous detection of excess at this wavelength is defined here as 
3cr above the expected white dwarf fiux. By demanding this level of excess, blackbody mod- 
els will suffice to predict the expected 4.5 /xm photospheric fiux of the white dwarf targets 
(Paper I). The absolute magnitude corresponding to the minimally reliable excess fiux from 
each white dwarf target, is then given by 

M4.5,„. = -2.51og(^(^) ) (2) 



where a is the total fiux error in Jy at 4.5 nm from Table HI d is the distance to the white 
dwarf in pc, and Fq = 179.7 Jy (ISpitzer Science Centerll2006a| ). Substellar cooling models 
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updated to include fluxes in the IRAC bandpasses were used to transfo rm the expected flu x 
into a mass for a given age (I. Baraffe 2007, private communication; iBaraffe et al.ll2003l ). 
Some representative, model predicted values for M4 at the ages of interest are given in 
Table [71 It is noteworthy that this analysis rules out unresolved sub-T dwarf companions at 
25 \ yhite dwarfs in Tab le M the known T dwarf sequence ends near M^,^^^ = 13.5 mag and 
T8 JPatten et al.lbooeh . and hmits at these stars reach 13.6 < M^,^^^ < 15.2 mag (typically 
M^,5^.m = 14.3 ±0.5 mag). 

Additionally, any spatially resolved point sources detected within several arcseconds 
of the white dwarf were photometrically examined for the possibility of companionship via 
their IRAC colors and any available ground-based photometric or astrometric data. Owing 
to the higher sensitivity of the 2 short IRAC wavelengths, some resolved substellar objects 
may not be detected at the 2 long IRAC wavelengths, resulting in a potential ambiguity for 
some visual companions. Generally speaking, no candidate companions were identified in 
this manner, but a few possibilties are discussed in §5.7. 



4.3. Resolved Companion Mass Limits from 7.9 yum Fluxes 



In order to detect T and sub-T dwarfs as spatially resolved companions, and to dif- 
ferentiate such objects from background point-like sources, they must be reliably detected 
at all 4 IRAC wavelengths. The 2 long wavelength IRAC channels in particular, together 
with the 2 short wavelength channels, provide unique information which should eliminate 
the color degeneracy between cool brown dwarfs and red extrag alactic (point-like) sources 
in the near-infrared and short wavelength IRAC channels alone (jPatten et al.ll2006l ). This 
fact appreciably limits any wide field IRAC searc h for T dwarf companions due to the lower 
sensitivity of the long wavelength IRAC channels (ISpitzer Science Centerll2006bl ). as detailed 
below. There are 7 white dwarfs from Table [3] which met the necessary criteria for such a 
search: 1) IRAC imaging of their surrounding fields in all 4 channels and; 2) a distance 
within approximately 20 pc. The entire known T dwarf sequence (down to T8) should be 
detected at all 4 wavelengths at these distances (jPatten et al.ll2006l ). 



Figure [9] shows the number of detected sources in each IRAC channel, as a function of 
magnitude, in the full IRAC fields of the 16 white dwarfs from Paper I which shared 600 
s integration times per filter, and identical 20-point dithering patterns. These sources were 
successfully detected and extracted photometrically by the IRAF tasks daofind and daophot. 
Based on the number of detections per magnitude bin (disregarding any trends in the number 
of source counts as a function of wavelength), it is clear that the 7.9 //m channel would limit 
any 4 channel IRAC survey for objects whose spectral energy distributions are not rising 
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towards longer wavelengths. 



4-3.1. Detectability of T and suh-T Dwarfs 



There exist a total of 58 white dwarfs which have been observed with IRAC at 7.9 
/im utilizing a common experimental design; 2 2 targets from the present w 'ork, 1 6 white 
dwarfs from Paper I, and 10 degenerates from lJura. Farihi. fc ZuckermanI (j2007al ). This 
white dwarf dataset allows an empirical assessment of the photometric sensitivity at this 
longest wavelength, and contains 35 unambiguous detections in that channel, with point 
sources as faint as 0.06 mJy reliably detected in all backgrounds. This find ing is consistent 
with: 1) the published sensitivities for IRAC (jSpitzer Science Centerll2006bl ): 2) calculations 
by the Sensitivity Performance Estimation Tool; and 3) the number of sources detected as 
a function of magnitude in Figure [H Therefore at my.g^m = 15.0 mag, or 0.064 mJy, point 
sources should be well detected regardless of background. 



^The T dwarf sequence down to spectral type T8 ends at My.g^m = 13.3 mag (iPatten et al 

20061 ). This corresponds to my.g^m = 14.8 mag at a distance of 20 pc and thus any widely 



separated T dwarf companions to the ci < 20 pc targets should be readily detected in this 
channel. Similar calculations in the other 3 IRAC channels estimate that a T8 dwarf should 
be correspo ndingly well detec ted out to; 20 pc at 5.7 yum, 60 pc at 3.6 /im, and 100 pc 
at 4.5 /im (jPatten et al.l 120061 ). For the white dwarf targets closer than 20 pc, the IRAC 
observations should be sensitive to small range of substellar companions with My.g^jm > 13.3 
mag, and potentially of a later spectral type than T. 

Owing to the nature of the dithering pattern, a further assessment must be made in 
regards to the effective field of view for the depth described above. For all but vMa 2, the 
medium sized cycling pattern was used, which should result in an effective coverage equal 
to the IRAC field of view minus about 25 pixels at each edge, or approximately 4.1' x 4.1', 
consistent with the analyzed images. For vMa 2, the effective field of view is about 12 pixels 
larger at each edge, or approximately 4.6' x 4.6', and the sensitivity was 0.75 mag less at 
each channel owing to a 150 s total integration time. 



4-3.2. Selection of T and suh-T Dwarfs 

For each white dwarf searched for wide T and sub-T dwarf companions, all 4 IRAC filter 
images were aligned and combined to create a single master coordinate image. The IRAF 
task daofind was executed on this master IRAC image to select point sources with counts at 
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or above 3a and the resulting coordinate list was then fed into daophot for each filter image 
in order to perform automated point spread function fitting photometry. 

Template point spread functions were created by running daophot on IRAC images of 
the Galactic component of the Spitzer Galactic First Look Survey to select approximately 
half to one dozen bright, unsaturated point sources in each filter. The magnitudes of these 
selected template soures were calculated by creating zero points which included the zero 
magnitude fluxes for IRAC, multiplication by the appropriate unit area on the array, and 
the necessary conversion of units. If the standard r = 10 pixel radius aperture is used for 
photometry, these zero points are (17.30, 16.81, 16.33, 15.69) mag at (3.6, 4.5, 5.7, 7.9) /xm. 

The extracted sources in each filter were cross-correlated using the master coordinate 
list, and 2 color-color diagrams were generated from the results; mse/^m — versus 
"^5.7Mm - mr.g^m, and mg.e^m - "^4.5Mm vcrsus m4.5^m " "^s-T/xm- In thcsc planes, T and 
sub-T dwarf candidates were sel ected by dernanding an object meet 3 criteria suggested 



by the IRAC colors presented in iPatten et al.l (120061 ): 0.2 < ms.g/xm — ^4.5Atm < 3.0, 0.0 < 



— fn7.9fira < 1-5, and —1.5 < m4.5^m — ?7^5.7/im < 1-0, with crror bars ignored in the 
first cut. All objects thus selected were examined individually; their images inspected and 
photometric data further evaluated. Although the mid-infrared colors of sub-T type objects 
are somewhat uncertain, model predictions yield colors which fit with the selection criteria 
above (D. Saumon 2006, private communication). Unfortunately, the earliest T dwarfs do 
not stand out strongly in I RAC color-color dia grams, and require near-infrared photometry 



to be clearly distinguished (IPatten et al.ll2006l ) 



No candidates which met all the criteria and passed critical examination were found 
in the IRAC fields of the 7 white dwarfs within d = 20 pc. While testing the color-color 
selection and extraction algorithm, the procedure was conducted at all 17 metal-rich white 
dwarfs in Table [21 A typical detection which met the color criteria had one or more of 
the following problems: 1) large photometric errors; 2) a location near the noisy edge of 
the image; 3) probable confusion with another source; 4) association with a known image 
artifact; 5) a substantially discrepant color-magnitude at the expected distance. Less often a 
detection would be a very red, unresolved extragalactic source whose nature was confirmed 
via existing optical astrometric and photometric catalogs. Table M lists the resulting upper 
limits on substellar mass companions achieved using this method, calculated by transforming 
^7.9fim = 15.0 mag to the expected Mr.g^m at the target distance, then using substellar 
cooling models to obtair i a mass from this fl ux at the appropriate total age (I. Baraffe 2007, 



private communication; iBaraffe et al.ll2003l ). 



To test the selection and extraction algorithm, the IRAC flelds of 3 T dwarfs were 
downloaded from the Spitzer archive (Program 35; PI Fazio) and the procedure was run 
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on images containing objects with spectral types T2.0, T5.0, and T8.0 (jPatten et al.ll2006l ). 
All 3 objects were selected correctly by the color-color cuts, and with modest photometric 
error bars indicating genuine detections. There is a non-zero probability that a bona fide 
cold brown dwarf companion escaped detection among our white dwarf target fields, despite 
the ability of daophot to spatially and photometrically deconvolve overlapping point sources. 
Based on the detection logs and the statistics from Figure [H in the two short wavelength 
IRAC filter images there were roughly 1000 sources per field with brightnesses greater than 
the completeness limit in those bandpasses. Taking a worst case scenario where all these 
sources represent potential spoilers yields a 2% probability of chance alignment within the 
relatively large, reduced image field of view. 



4.3.3. The Hyades 

Figure dU] plots color-magnitude diagrams for all detected point-like sources in the fields 
of the 6 Hyades white dwarfs observed in both of the short w avelength IRAC ch annels. 
Included in the plot is the expected T dwarf sequence (T1-T8; IPatten et allbood ) at the 



46 pc distance to the open cluster. The cooler part of the T dwarf sequence appears to 
standout from most field objects. All sources with ms.g/im — > 1-0 were investigated 

individually, revealing a few extragalactic sources, spurious detections near bright stars or the 
edge of the mosaic, and sources with large photometric errors. No reliable candidates near 
the T dwarf sequence were identified. Since T and sub-T dwarfs should not be detected at 
the distance to the Hyades in the 2 long wavelength IRAC channels, this limits what can be 
done with single epoch IRAC data to constrain wide, methane-bearing su bstellar companions . 



Moreover, because very cool brown dwarfs will have K — m^ ^^^ > 1 (jPatten et al.ll2006l ). 
they should not be detected in 2MASS, as any such companions in the Hyades would have 
> 16 mag. Therefore, further analysis of these datasets can only be be achieved with 
deep near-infrared imaging from the ground or proper motion analysis with a second epoch 
IRAC observation. 



5. LIMITS ON MASSIVE PLANETS AND COLD BROWN DWARF 

COMPANIONS 

Clearly, none of the 32 stars in Figures [TH8] display reliably measured photometric excess 
at 4.5 yum. When combined with similar results for the white dwarfs analyzed in Paper 
I, there are a total of 47 observed degenerates which reveal no evidence for unresolved, 
cold brown dwarf or massive planetary companions. This is a striking result, especially 
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considering the large number of relatively young white dwarfs where massive planets would 
be detectable. There are 34 white dwarfs for which the IRAC 4.5 fim observations were 
sensitive to unresolved planetary mass companions in the range 2 — 13 Mj, and 10 white 
dwarfs for which the data were sensitive to brown dwarf companions in the range 14 — 20 Mj, 
according to subst ellar cooling models used here (I. Baraffe 2007, private communication; 



Baraffe et al.l 120031 ). 



5.1. The Influence of Stellar Evolution on Massive Planets 



The aperture photometry and image analysis places limits on the presence of unresolved 
or partially resolved companions out to approximately 5 pixels or 6", the largest aperture 
used for photometry (Paper I). This angle on the sky corresponds to a several hundred AU 
region around the target white dwarfs, which lie at typical distances in the range d ~ 20 — 100 
pc. However, any planets which formed in the 5 — 50 AU range during the main sequence 
phase, should now be locate d further out due to orbital expansion via mass loss during 
the asymptotic giant branch (Farihi^ Hoard, fc Wachtei [2OO6 : Farihi. Zuckerman. fc Becklin 



2005bl : iBurleigh. Clarke, fc Hodgkin 



2002 



Zuckerman fc Becklin 



1987 



Jeanslll924j ). yet still 



within a region to which the IRAC photometry is sensitive. While a typical expansion factor 
(equal to the ratio of the main sequence progenitor mass to the white dwarf mass) is 2 — 3 for 
F-type stars, the massive young degenerates studied here are likely to produce much larger 
increases, up to a factor of 6 or so. 

Tidal interactions should be relatively strong for massive planets orbiting intermediate 
mass stars, whose radii can become as large as several AU at the tip of the asymptotic giant 
bran ch. For example, a 5 M p^ main sequence star should grow to a maximum radius near 5 
A U (IVillaver fc Livid 120071 ). directly engulfing any planets in that range. Using Equation 6 
of iRasio et al.l (119961 ) for such a star, and assuming a convective envelope mass near 3 Mq, a 
10 Mj planet at 10 AU should tidally d ecay within the 1 Myr lifetime of its asymptotic giant 
branch host (IVassiliadis fc Woodlll993l ). Although this calculation is likely oversimplified, it 
is instructive; eschewing direct engulfment is not sufficient for a planet to survive the post- 
main sequence. Rather, the more massive the planet, the more it will induce tides in the 



giant star and the more likely it will experience orbital decay and be destroyed (IRasio et al. 



19961 ). One competing factor is the fact that by the time an asymptotic giant has reached 
it maximum radius there has been significant mass lost and the orbits of any planets should 
already have expanded commensurately. However, a complete, time-dependent treatment of 
these competing forces - orbital expansion due to mass loss versus orbital decay due to tidal 
forces - on planets during the asymptotic giant branch has not been carried out. Such a 
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study would be highly valuable but is beyond the scope of this paper. 



5.2. Massive Planets and Brown Dwarfs at Evolved Stars 



Recent work indicates a substantial percentage of substellar, radial velocity companions 
to first ascent giant stars are potentially or certainly above the deuterium burning minimum 
mass. Presentl y 4 of 13 or 31% of known substellar companions to giant star s have masses 



abov e 13 Mj flLiu et all 120071 : iLovis fc Mavorl 120071 : iNiedzielski et all 120071 : iHatzes et al. 



20051 ) . Two items of interest for white dwarf planetary system studies emerge from these 
results at giant stars. First, the Msini distribution of close substellar companions to inter- 
mediate mass, evolved stars is markedly different than for solar-type main sequence stars. 
Second is the fact that all 4 brown dwarf hosting giants have 2 Mq < M < 4 Mq, and are 
thus related to the population of white dwarfs studied in this paper. 



With orbital semim a.ior axes 0.5 AU < a < 2.5 AU (ILiu et al.ll2007l : lLovis fc Mavorll2007l : 
Niedzielski et al.ll2007l ). all of the known substellar radial velocity companions to first ascent 
giant stars risk destruction during the ensuing asymptotic giai it phase. However, the most ap- 
paren tly brown dwarf-like companions have M sin i ^ 20 Mj (ILiu et al.ll2007l : iLovis &: Mayor 
20071 ): potentially massive enough to eject the giant stellar envelope as have the 50 — 60 M- 



substellar companions to the white dwarfs WD 0137—049 and GD1400 (IBurleig. 



h et al. 2006 



Maxted et al.ll2006l : iFarihi. Zuckerman. fc BecklinI l2005bl : iFarihi fc Christophei] 120041 11 The 
close, P 2 hr substellar companion to WD 0137—049 is thought to have ejected the dense 
first ascent giant envelope of its host, as evidence by the fact the white dwarf is a low mass, 
helium core degenerate. This is not the case for GD 1400B, whose degenerate host is a 
typical, carbon-oxygen core white dwarf, but its P ^ 10 hr period (Burleigh et al. 2008, in 
preparation) indicates likely prior orbital decay due to the ejection of the asymptotic giant 
envelope. By inference, both these substellar survivors would have originally orbited within 
roughly 1 AU of their host intermediate mass stellar progenitors. 



5.3. Formation, Persistence, and Sacrifice 



The results of this IRAC white dwarf study may imply that close massive planetary and 
brown dwarf companions to intermediate mass stars do n ot typically surv i ve th e asymptotic 
giant branch. Robust statistics are not yet available, but iLovis &: Mayorl (120071 ) estimate at 
least 3% of stars with M > 1.8 Mq mass stars host Msini > 5 Mj companions, and as 
discussed in §5.2, about 1/3 of these are brown dwarfs (according to the lAU definition). 
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Hence, a reasonably optimistic expectation for an IRAC search of 46 white dwarfs would 
be 1 or 2 detections. This negative result suggests the possibility that a higher (substel- 
lar) companion mass is required to survive the entire post-main sequence - specifically the 
asymptotic giant phase - at orbital separations of a few to several AU. 

However, recent evidence has suggested a mechanism which may allow planets to survive 
within the inner regions through both giant phases; sacrifice of the innermost components. 
Analogous to the formation of several types of post-main sequence bina ries, substellar com- 
panions cari dyna mically eject the (first or second phase) giant envelope (INelemans &: Tauris 
19981 : ISokerlll998l ). thereby shielding any remaining components in a planetary system. This 
process - common envelope evolution - is responsible for the close orbits of numerous low 



mass stellar and substellar companions to white dwarfs (IFarihi. Hoard, fc Wachterl 12006 



Schreiber &: Gansickdl2003l ). Essentially, unstable mass transfer from the giant results in a 
frictional exchange of orbital energy between the secondary and the envelop e material, result- 



ing in an orbit decrease and an efficient envelope ejection (iPaczynskil 
Mj, a = 1.7 AU planet recently detected at the sdB star V391 Pegasi (ISilvotti et al.ll2007l ) has 



1976). The Msin; 



survived a first asceiit giant phase involving significa r it mass loss, possibly due to such an in- 
teraction (IHan et al.l 120021 : iNelemans fc Taurislll998l : ISoken 119981 ). Sim ilarly, the Msini ^ 2 
Mj, a = 2.5 AU planet candidate at the pulsating white dwarf GD 66 (jMuUally et al.ll2008l ). 
may owe its inner region survival to the sacrfice of closer planets. 

While intermediate separation planets may have indeterminate fates, any massive plan- 
ets originally orbiting at a > 5 AU, should now be located at tens to hundreds of AU in 
the white dwarf phase, havi ng essentially eluded a ny substantial post-main sequence inter- 
action with their host star (jVillaver fc Livid 120071 ). The results of this IRAC search imply 
that massive, M > 10 Mj planets and brown dwarfs form rarely (/ < 3%) at these wide 
separations. 



5.4. Planets in the Hyades 



The massive planet recently found at the Hyades giant e Tauri (ISato et al.ll2007l ) is very 
interesting in light of the 7 classical white dwarf Hyads surveyed with IRAC - would such an 
planet have been detected if it persisted into the white dwarf evolutionary phase? The most 
likely mass of the planet at e Tauri is near 10 Mj, and the IRAC 4.5 photometry was 
sensitive to objects of this mass at virtually every Hyades target, although in reality such 
a detection might prove more or less difficult due to differences from the model predictions 
used here. However, it is uncertain whether this Hyades planet at a = 1.9 AU will survive 
the asymptotic giant phase of its host, since the maximum radius of the star will reach at 
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least 3 AU fIVillaver fc Livioll2007h . 

A clear advantage of the IRAC search of the Hyades white dwarfs is insensivitiy to 
orbital separation or inclination, parameter spaces which limit radial velocity and direct 
imaging searches for planets. The resulting substellar mass sensitivities achieved here for 
the Hyades are cor nparable to those pro duced via direct imaging with HST / NICMOS - 
also around 10 Mj (IFriedrich et al.ll2005l ) - but for any orbits out to a ~ 250 AU. Only very 
widely separated massive planets should have escaped detection. 



5.5. Cold Brown Dwarfs 



The previous sections have covered discussions of mass, but not of temperature. For 
27 targets in the survey, the IRAC 4.5 /im photometry was sensitive to the entire known T 
dwarf sequence, independent of the corr esponding masses, and in a few cases more than a 
full magnitude fainter than a T8 dwarf (jPatten et al.ll2006l ). In fact, with the exception of 
a single target (G21-16; see §5.7) whose IRAC image was confused with other sources, this 
survey rules out brown dwarf companions down to 25 Mj within a few hundred AU of all 
white dwarf targets, implying a brown dwarf companion fraction less than around 2%. If the 
minimal sensitivity achieved here f or these 46 targets i s matched by similar IRAC 4.5 /im 
results for 124 nearby white dwarfs (iMuUally et al.ll2007l ). this fraction could be smaller than 
0.6%, and consistent with previous white dwarf studies whi c h were sensitive to somewhat 
higher temperatures and masses (L dwarfs: iHoard et al.ll2007l: iFarihi. Becklin. fc Zuckerman 

inlll987h. 



2005a; Farihi 2004: Wachter et al. 2003: Zuckerman k Beckli: 



5.6. No Evidence for Merger Products 



Although not the focus of this study, the results give no indication of lasting merger 
products at any of the 10 addit ional massive white dwarfs observed speci fically for this study. 
Combined with the results of [Hansen. Kulkarni. fc Wiktorowiczl (120061 ). there is as yet no 
indication of disks (and subsequent pollution) or reforged planets at nearly 2 dozen massive 
white dwarfs which could conceivabl y have formed as a merger of lower mass degenerates 
( iLiebert. Bergeron. &: Holberg||2005al ). Given the recent results on externally-polluted white 
dwarfs with debris disks, it is possible that any dust disks formed via white dwarf mergers 
would dissipate rapidly, the particles becoming gaseous through mutual collisions within the 
disk (Paper I). 



- 16 - 



5.7. Individual Objects 

0046+051 (vMa 2) At 4.4 pc, this degenerate represents a unique and advantageous 
hunting ground for planets and planetary system remnants. This cool, helium atmosphere, 
metal-rich white dwarf has been externally polluted by interstellar or circumstellar mat- 
ter. Previous ground- and space-based mid-infrared ima ging and photo metry have ruled 
out the presence of a substellar companion suggested bv iMakarovl (l2004r). down to T^ff > 



500 K and corresponding to aroun d 25 Mj at 5 Gyr (IFarihi. Becklin. fc Macintosh! 12004 



Burrows. Sudarsky. fc Lunind 120031 ). The IRAC 4.5 jj,m photometry of vMa 2 and the mod- 
els used here rule out the presence of a companion as cool as 400 K, and a mass close to the 
deuterium burning limit at 4.4 Gyr. Furthermore, the IRAC 4 channel color-color search for 
resolved substellar companions rules out the presence of any widely bound object as cool as 
Teff ~ 550 K within r 1200 AU of vMa 2. Models predict this should correspond to a mass 
of 25 Mj at the age of this well studied white dwarf. 

Deep ground-based J-band imaging observations have ruled out widely-bound planetary 
mass companions to vMa 2 as small as 7 Mj at orbital separations near 10 — 200 AU 



(IBurleigh et al.ll2008l ). There exist 2 epochs of IRAC 4.5 /xm imaging of vMa 2 in the Spitzer 
archive, separated by 2.1 yr and clearly revealing 6.2" of proper motion upon blinking the 
aligned frames. There are no field objects comoving with the white dwarf, ruling out well 
detected, resolved objects with m4.5^jn = 16.7 mag (see Figure [H]) or 36 /xJy as companions 
within 1200 AU of vMa 2. At the 4.4 pc distance and 4.4 Gyr age of vMa2, models predict 
that this very sensitive observational limit corresponds to a mass of 4 Mj and a temperature 
of just Tefj ~ 200 K; signficantly lower than the known T dwarfs and only 40 K warmer than 
Jupiter(!). 

0236+498 (EUVE J0239-F50.0) Very little reliable photometry exists on this hot and 
faint degenerate, making it difficult to properly calibrate its spectral energy distribution (see 
Figure [2]). The 2MASS J-band flux is likely the most reliable photometric data available, 
while the iJ-band flux has a large associated error. The white dwarf is clearly detected in the 
short wavelength IRAC images, but significant uncertainty exists at the long wavelengths, 
where the fiux in the photometric aperture could be due to a background source or noise. 
Unfortunately, complications plague the short wavelength IRAC observations; a column pull- 
down artifact at the position of the white dwarf and a photometrically overlapping point 
source located 3.6" away. The IRAF routine daophot was used to deconvolve the fiux of the 
white dwarf and the nearby point source, after manually correcting the column pull-down 
artifact by adding the median value of nearby columns. Despite these efforts, it is possible 
or even likely the photometry of the white dwarf is contaminated. If the 2MASS if-band 
flux and IRAC photometry are somewhat accurate, then the white dwarf may possess an 
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excess. Ground-based near-infrared photometry is needed for this source. 



0325—857AB (LB 98 02AB) This visual pair has long been a binary suspect (IHansen. Kulkarni. &: Wikt 



20061 : iBarstow et al.lll995l ) but its physical companionship is confirmed here for the first time. 
Previous arguments for binarity relied on space density, which is insufficient to firmly es- 
tablish companionship, especially in a rare system such as this, where the more massive of 
the pair white dwarfs is the hotter and apparently younger star. Common proper motion is 
demonstrated by measuring the astrometric offsets of the white dwarfs and 20 background 
point sources between 2 epochs of SuperCOSMOS data. Two scans of UKST photographic 
plates were used for this purpose, one image taken in 1976.7 and another from 2001.9, pro- 
viding a 25.3 yr baseline. Using the IRAF routine geomap to measure both the motion of 
the white dwarfs and the random centroiding errors from the positions of the background 
sources, the component proper motions are: /i^ = 0.074" yr^^ at 6a = 244° and /i^ = 0.083" 
yr~^ at 9b = 242°, with = 0.012" yr~^ and cxe = 6°. The separation between primary and 
secondary from the 2MASS J-band image is 6.9" at position angle of 325°. 

1455+298 (G166-58) This metal-rich white dwarf displays excess flux in the 2 long 
wavelength IRAC channels, but not at 4.5 /im (Paper I). Therefore, it was included in the 
photometric analysis for unresolved substellar companions. 

2326+049 (G29-38) Not included in the analysis for unresolved companions due to its 
large photometric excess at all IRAC wavelengths (Paper I). However, the field surround- 
ing G29-38 was searched for widely-bound, cold substellar companions using the method 
described in §4.3, revealing no candidates down to 25 Mj within 3500 AU. 



6. CONCLUSION 

A relatively comprehensive look at both young and older white dwarfs with IRAC re- 
veals no promising evidence for massive planets or cold brown dwarf companions at orbital 
separations within a few hundred AU. By conducting a search for substellar companions via 
excess emission in the mid-infrared, the data cover the widest range of orbital phase space 
possible; nearly all possible semimajor axes, eccentricities and inclinations, including the 
missed middle ground (5 — 50 AU) between radial velocity and direct imaging searches. The 
uncovered range of possible orbits lies roughly beyond a few hundred AU. 

While it is somewhat uncertain how substellar objects evolve (dynamically or otherwise) 
within 5 — 10 AU of mass losing asymptotic giant stars, avoidance of direct engulfment may 
not suffice for ultimate survival to the white dwarf phase. However, massive inner planets 
could act as sacrificial guardians for remaining outer planetary system components. Any 
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planets and brown dwarfs outside of this region should be relatively unaffected, and the 
IRAC results demonstrate a dearth of cold substellar companions down into sub-T dwarf 
temperatures. The 2 — 3 Mj planets suspected at the sdB star V 391 Pegasi and the white 
dwarf GD 66 are just beyond the reach of this survey. 

The negative results for the Hyades and comparably aged white dwarfs yield upper mass 
limits at or somewhat below 10 Mj for objects which may have formed around 3 Mq main 
sequence stars. Similar limits at the Pleiades and analogously hot and massive field white 
dwarfs provide the first evidence that massive planets are not commonly formed or do not 
survive the post-main sequence evolution of intermediate mass stars with M > 4 M©. This 
latter result was only achievable via observations of white dwarfs, as their main sequence, 
B-type progenitors are not amenable to other planet detection techniques. 

The lack of 4.5 /x m excess at all white dwarf targets, especially when combined with 



similar IRAC searches (jMuUally et al.ll2007l ). confirms that L and T-type brown dwarf com- 
panions are rare (/ < 0.6%) within a few hundred AU, down to masses near the deuterium 
burning limit. These results suggest the possibility that the lowest mass companions, and 
especially planets, orbiting intermediate mass stars may be altered or destroyed prior to or 
during the post-main sequence, or are (more likely) too cold to directly detect with current 
facilities. 
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Table 1. Young Cluster White Dwarf Targets 



WD# 


Name 


Teff (K) 


V (mag) 


M (Mq) 


Referen 


0349+247 


EG 25 


31700 


16.64 


1.09 


1,2,3 


0352+096 


EG 26 


14800 


14.52 


0.71 


1,2,4 


0406+169 


EG 29 


15200 


15.35 


0.80 


1,2,4 


0415+271^ 


EG 265 


11500 


15.00 


0.55 


5,6,7 


0421+162 


EG 36 


19600 


14.29 


0.68 


1,2,4 


0425+168 


EG 37 


24400 


14.02 


0.70 


1,2,4 


0431+126 


EG 39 


21300 


14.24 


0.65 


1,2,3 


0437+138 


GR 316 


15300 


14.93 


0.68 


1,2,7 


0438+108 


EG 42 


27400 


13.86 


0.75 


1,2,3 



'^EG 265 i s a member of the Hyades cluster or supercluster ( ReidI 
Eggeij|l984h . 



1993 



(7) lUpgren. Weis. fc Hanson 



References. — (1) 


Claver et al. 2001; (21 


Bergeron. Liebert. & FuUbright 


1995a; (3) Cheselka et al. 


1993 


C4) Eggen & Greenstein 


1965; (5) Bergeron et al. 2004; (6) Reid 


1993 



198J 
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Table 2. Young Field White Dwarf Targets 



WD# 


Name 


Teff (K) 


V (mag) 


M (Mo) 


References 


0001+433 


EUVE 


42400 


16.8 


1.37 


1,2,3 


0136+251 


PG 


39400 


15.87 


1.32 


1,3,4 


0235-125 


PHL 1400 


32400 


14.98 


1.03 


1,2,5 


0236+498 


EUVE 


33800 


15.8 


0.94 


6,7 


0325-857A 


LB 9802A 


16200 


14.11 


0.85 


8,9 


0325-857B 


LB 9802B 


33800 


14.90 


1.33 


1,8,10 


0346-011 


GD 50 


43200 


14.04 


1.37 


1,2 


0440-038 


EUVE 


65100 


16.7 


1.33 


1,2,7 


0518-105 


EUVE 


33000 


15.82 


1.07 


1,2,3 


0531—022 


EUVE 


29700 


16.20 


0.97 


6,7 


0652-563 


EUVE 


35200 


16.6 


1.18 


1,2 


0730+487 


GD 86 


15500 


14.96 


0.90 


12 


0821-252 


EUVE 


43200 


16.4 


1.21 


13 


0914-195 


EUVE 


56400 


17.4 


1.33 


1,2 


1022-301 


EUVE 


35700 


15.9 


1.27 


1,6,7 


1440+753^ 


EUVE 


35000 


15.22 


1.06 


3,13,14 


1529-772 


EUVE 


51600 


16.9 


1.24 


2,13 


1543-366 


EUVE 


45200 


15.81 


1.17 


6,15 


1609+631 


PG 


30400 


16.68 


1.05 


4,6 


1642+413^ 


PG 


26500 


16.21 


0.96 


4,6 


1658+440 


PG 


30500 


15.02 


1.41 


1,4 


1740-706 


EUVE 


46800 


16.51 


1.18 


1,13,15 



''1440+753 is a close double degenerate 



'iFinlev. Koester. fc Basril (Il997l) lis t 0.96 M© (109pc) for 1642+413, while 
Liebert. Bergeron, fc Holber j |2005al ) give 0.79 Mq (145 pc) 



References. — (1) Hansen Kulkarni. & Wiktorowicz 20061 (21 


Vennes et al. 1997; (3) Marsh et al. 1997; (4) Liebert. Bergeron. & Holberg 


2005a; (5) I 


DuDuis. Vennes. & Chaver 2002; (6) 


Finlev. Koester. & Basri 


1997 


; (7) McCook & Sion 2006; (8) Barstow et al. 


1995; (9) Ferrario et al. 


1997 


; (10) 


Vennes et al. 


20031: (12) iBergeron. Saffer. & LiebertJ 19921 


(13) Vennes 


1999; (U) 


Vennes. Ferrario. & Wickramasinghel 1999; fl5) 


Vennes et al. 


1996 



Table 3. Metal-Rich Field White Dwarf Targets 



WD# 


Name 


Toff (K) 


V (mag) 


M (Mq) 


References 


0032-175 


G266-135 


9240 


14.94 


0.60 


1,2 


0046+051 


vMa 2 


6770 


12.39 


0.83 


3 


0235+064 


PG 


15000 


15.5 


0.61 


4 


0322-019 


G77-50 


5220 


16.12 


0.61 


5,6 


0846+346 


GD 96 


7370 


15.71 


0.59 


1,7 


1102-183 


EC 


8060 


15.99 


0.60 


1,7 


1124-293 


EC 


9680 


15.02 


0.63 


5,8 


1204-136 


EC 


11500 


15.67 


0.60 


1,9 


1208+576 


G197-47 


5880 


15.78 


0.56 


3 


1344+106 


G63-54 


7110 


15.12 


0.65 


3 


1407+425 


PG 


10010 


15.03 


0.73 


10 


1455+298 


G167-8 


7390 


15.60 


0.58 


3,10 


1632+177 


PG 


10100 


13.05 


0.58 


10 


1633+433 


G180-63 


6690 


14.84 


0.72 


3,10 


1826-045 


G21-16 


9480 


14.58 


0.57 


3 


1858+393 


G205-52 


9470 


15.63 


0.60 


1,7 


2326+049 


G29-38 


11600 


13.04 


0.69 


7,11 



Note. — For those stars with no spectroscopic or trigonometric mass- 
radius estimate, log g = 8.0 was assumed. 



References. — (1) Zuckerman et al. 


2003; (2) 


MermiUiod 


1986: (3) Bereeron. Leeeett. & RuizI 2001; (4) This work: (5) 


Bergeron. Ruiz. & Leesett 1997 


; (6) 




Smart et al. 2003: (7) 


McCook & Sion hood (8) 


Koester et al. 




2001: (9) 


Salim & Gould 


2003: (10) Liebert. Bergeron. & Holbere 


2005a; fll) 


Bergeron et al. 


1995d| 



-28- 



Table 4. IRAC Fluxes for White Dwarf Targets 



vv -L'tr 


fi..™ { //.7vi 


Wa k..™ (//.7vi 


T^K. '7,,™ ( /y„Tvl 


p7 o,<™ f//,.7v) 


J- XL/^-'XXXXV^ 


0001+433 


18 + 6 


13 + 7 


30"* 


29^^ 


14.0 


0046+051 


8040 ± 400 


5360 ± 270 


3680 ± 190 


2080+ 110 


14.0 


0136+251 


46 ±6 


26 + 7 


34^ 


35^^ 


14.0 


0235-125 


111 + 8 


63 + 7 


36^ 


35^^ 


14.0 


0236+498 


107+11 


72 + 9 


20^ 


27^^ 


11.0 


0325-857A 


396 ± 20 


238 ± 14 


139 ± 27 


87 + 25 


14.0 


0325-857B 


149 + 9 


92 + 8 


53 + 28 


33^^ 


14.0 


0346-011 


263 ± 14 


164+11 


145 ± 34 


47a 


14.0 


0349+247 


27 + 3 


15 + 4 


22"* 


39*^ 


11.4 


0352+096 


295 ± 15 


168 + 9 


99 + 21 


44 + 28 


11.4 


0406+169 


134 + 7 


69 + 5 


55 + 20 


61 + 34 


11.4 


0415+271 


245 ± 13 


148 + 8 


100 + 19 


68 + 23 


11.4 


0421+162 


275 ± 14 


163 + 9 


121 ± 20 


58 + 23 


11.4 


0425+168 


313 + 16 


185 + 10 


133 ± 21 


101 + 27 


11.4 


0431+126 




179+10 




5P 


14.0 


0437+138 


186 + 10 


107 + 6 


62 + 17 


42^ 


11.4 


0438+108 




211 + 12 




76 + 38 


14.0 


0440-038 


20 + 5 


19 + 6 


23^ 


2P 


14.0 


0518-105 


51 + 5 


33 + 7 


32^ 


31^^ 


14.0 


0531-022 


47 + 4 


26 + 5 


20*^ 


37^^ 


11.4 


0652-563 


16+11 


18 + 9 


32"* 


30*^ 


14.0 


0730+487 


218+11 


127+7 


98 + 22 


66 + 24 


10.5 


0821-252 


26 + 10 


13 + 7 


21"* 


20^^ 


11.0 


0914-195 


13 + 5 


9 + 6 


33"* 


32*^ 


14.0 


1022-301 


33 + 6 


22 + 7 


30^ 


31*^ 


14.0 


1440+753 


88 + 6 


58 + 5 


27+16 


16^^ 


11.0 


1529-772 


22 + 5 


12 + 4 


20^ 


18^ 


12.4 


1543-366 


51 + 5 


31 + 5 


23"* 


23*^ 


11.4 


1609+631 


24 + 3 


13 + 3 


15*^ 


22=^ 


11.0 


1642+413 


38 + 3 


21 + 4 


16'^ 




11.4 


1658+440 


131 + 8 


81 + 7 


56 + 27 


28^ 


14.0 


1740-706 


25 + 4 


14 + 4 


16"* 


19*^ 


14.0 



Upper limit. 



Note. — Photometric errors and upper limits are described in §3.1. 
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Table 5. Target Ages and Upper Mass Limits for Unresolved Companions 



WD# 


tms (Gyr) 


^wd (Gyr) 


(Gyr) 


d (pc) 


Mi.si^m (mag) 


Mass (Mj) 








Young 


Cluster White Dwarfs 






0349+247 






0.1 


132 


12.3 


9 


0352+096 






0.6 


46 


13.7 


10 


0406+169 






0.6 


46 


14.4 


7 


0415+271 






0.6 


46 


13.9 


9 


0421+162 






0.6 


46 


13.7 


10 


0425+168 






0.6 


46 


13.6 


10 


0431+126 






0.6 


46 


13.6 


10 


0437+138 






0.6 


46 


14.2 


8 


0438+108 






0.6 


46 


13.4 


11 






Youn; 


I Field White Dwarfs 






0001+433 






0.1 


96 


12.4 


9 


0136+251 






0.1 


73 


13.0 


6 


0235-125 






0.1 


65 


13.3 


5 


0236+498 


0.1 


0.1 


0.2 


107 


11.9 


11 


0325-857A 


0.1 


0.3 


0.4 


35 


13.9 


6 


0325-857B 






0.4 


35 


14.5 


5 


0346-011 






0.1 


30 


14.5 


3 


0440-038 






0.1 


134 


11.9 


10 


0518-105 






0.1 


92 


12.5 


8 


0531-022 


0.1 


0.1 


0.2 


107 


12.6 


9 


0652-563 






0.1 


119 


11.7 


11 


0730+487 


0.1 


0.4 


0.5 


39 


14.5 


6 


0821-252 






0.1 


105 


12.2 


10 


0914-195 






0.1 


175 


11.3 


12 


1022-301 






0.1 


61 


13.4 


5 


1440+753 






0.1 


101 


12.7 


7 


1529-772 






0.1 


137 


12.3 


9 


1543-366 






0.1 


111 


12.5 


8 


1609+631 






0.1 


134 


12.6 


8 


1642+413 


0.1 


0.1 


0.2 


109 


12.8 


8 


1658+440 






0.1 


32 


14.8 


2 
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Table 5 — Continued 



WD# 


tms (Gyr) 


twd (Gyr) 


(Gyr) 


d (pc) M4. 


5^m (mag) 


Mass (Mj) 


1740-706 






0.1 


76 


13.5 


5 


Metal-Rich Field White Dwarfs 


0032-175 


1.5 


0.7 


2.2 


31 


13.9 


18 


0046+051 


0.2 


3.7 


3.9 


4.4 


15.2 


13 


0235+064 


1.2 


0.2 


1.4 


70 


13.3 


20 


0322-019 


1.2 


4.4 


5.6 


17 


15.1 


17 


0846+346 


1.9 


1.4 


3.3 


30 


14.6 


17 


1102-183 


1.5 


1.1 


2.6 


40 


14.2 


17 


1124-293 


0.9 


0.6 


1.5 


34 


14.3 


12 


1204-136 


1.5 


0.4 


1.9 


62 


13.5 


20 


1208+576 


4.0 


2.5 


6.5 


20 


14.7 


25 


1344+106 


0.7 


1.9 


2.6 


20 


14.7 


14 


1407+425 


0.3 


0.9 


1.2 


33 


14.6 


10 


1455+298 


2.4 


1.4 


3.8 


36 


14.1 


25 


1632+177 


2.4 


0.6 


3.0 


16 


14.1 


20 


1633+433 


0.3 


2.7 


3.0 


15 


14.8 


14 


1826-045 


3.0 


0.8 


3.8 


29 


13.0 


35 


1858+393 


1.5 


0.7 


2.2 


45 


14.2 


16 



'^The total age of the Pleiad EG 25 (Table [T]) and the similarly hot and massive field white 
dwarfs (Table[2]) is taken to be 0.125 Gyr (see §4.1). The seventh column lists the companion 
mass upper limit. 
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Table 6. Target Ages and Upper Mass Limits for Resolved Companions 



WD# 


tnis (Gyr) 


twd (Gyr) 


T (Gyr) 


d (pc) 


Mrsfim (mag) 


Mass (Mj) 


0046+051 


0.2 


3.7 


3.9 


4.4 


16.0 


25 


0322-019 


1.2 


4.4 


5.6 


17 


13.8 


65 


1208+576 


4.0 


2.5 


6.5 


20 


13.5 


70 


1344+106 


0.7 


1.9 


2.6 


20 


13.5 


60 


1632+177 


2.4 


0.6 


3.0 


16 


14.0 


50 


1633+433 


0.3 


2.7 


3.0 


15 


14.1 


50 


2326+049 


0.4 


0.5 


0.9 


14 


14.3 


25 



Note. — The exposure time for 0046+051 was 150 s versus 600 s for the other targets, 
resulting in an overall sensitivity about half of that calculated in §4.3.1. 



Table 7. Representative 4.5 /im Absolute Magnitudes for Substellar Objects 



Mass (Mj) 


0.1 Gyr 


0.5 Gyr 


1.0 Gyr 


5.0 Gyr 


2 


14.8 


16.7 


17.7 


21.1 


5 


13.2 


14.9 


15.7 


17.9 


10 


12.0 


13.5 


14.3 


16.1 


15 


11.1 


12.7 


13.5 


15.3 


20 


10.9 


12.2 


12.9 


14.7 



Note. 



Entries are in magnitudes 
The table is based 



(jSpitzer Science Centeij l2006al ). 
on the models of I. BarafFe (2007, private communication 



Baraffe et all 120031 ) 
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1.000 F 



EUVE J0003%3.5 



100.0 F 



T = 42400 K 



0.010 : 



1.0 10.0 

X {fj.m) 




10.00 



:-PG. 0136+251 



0.01 

0, 



10.00 



T . = 39400 K 



: PHL 1400 



I. T I 



1.0 10.0 



^ 0.01 
100.0 



T = 32400 K 



I -.J 



1.0 10.0 
X (/im) 



Fig. 1.— Spectral energy distribution of EUVE J0003+43.5, vMa 2, PG 0136+251, and 
PHL 1400. Downward arrows represent upper limits (§3.1). 
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10.00 



0.10 



EUVE J0239+50.0 



10.00 : 



33800 K 

IRAC BEAM CONTAMINATED 



0.10 - 



0.01 



LB 9802A 



T = 16200 K 



1.0 10.0 



100.0 0.1 



1.0 10.0 100.0 



10.00 



: LB 9802B ' 



10.00 



T . = 33800 K 



: GD 50 



0.01 L_ 
0.1 



1.0 10.0 



^ 0.01 L_ 

100.0 0.1 



T = 43200 K 



1.0 10.0 



Fig. 2.— Spectral energy distribution of EUVE J0239+50.0, LB 9802A, LB 9802B, and GD 
50. Downward arrows represent upper limits (§3.1). 
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10.00 F 



: EG 25 



10.00 F 



T = 31700 K 



0.10 : 



0.01 L_ 
0.1 



V. 



1.0 10.0 

X {fj.m) 



0.10 - 



^ 0.01 
100.0 



EG 26. 



T = 14800 K 



1.0 10.0 



10.00 F 



: EG 29 



10.00 



T . = 16200 K 



: EG 265 



0.01 L 

0.1 



T = 11500 K 



1.0 10.0 



^ 0.01 L_ 

100.0 0.1 



1.0 10.0 100.0 



Fig. 3. — Spectral energy distribution of EG 25, EG 26, EG 29, and EG 265. Downward 
arrows represent upper limits (§3.1). 
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10.00 F 



: EG 36 



10.00 F 



T = 19600 K 



0.10 : 



1.0 10.0 

X {fj.m) 



0.10 - 



0.01 L 
0.1 



EG 37 



T = 24400 K 



1.0 10.0 



10.00 F 



: EG 39 



10.00 



T . = 21300 K 



: GR 316 



I: 



0.01 L 

0.1 



T = 17000 K 



1.0 10.0 



^ 0.01 L_ 

100.0 0.1 



1.0 10.0 100.0 



Fig. 4. — Spectral energy distribution of EG 35, EG 37, EG 39, and GR 316. Downward 
arrows represent upper limits (§3.1). 
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10.00 F 



: EG 42 



10.00 F 



T = 27400 K 



0.10 : 



1.0 10.0 

X {fj.m) 



0.10 - 



0.01 L 
0.1 



■EUVE J0443-03.7 



T = 66100 K 



1.0 10.0 
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: EUVE J0521-10.4 



10.00 
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0.01 L_ 
0.1 
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1.0 10.0 



^ 0.01 L_ 

100.0 0.1 
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Fig. 5.— Spectral energy distribution of EG 42, EUVE J0443-03.7, EUVE J0521-10.4, 
and EUVE J0534— 02.2. Downward arrows represent upper limits (§3.1). 
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1.000 F 



: EUVE J0653^6.4 



10.00 F 



T = 36200 K 



0.010 : 



I I 



1.0 10.0 

X {fj.m) 




1.000 



: EUVE J0823-25.4 



1.000 



T . = 43200 K 
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^ 0.001 L_ 
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1.0 10.0 100.0 
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Fig. 6.— Spectral energy distribution of EUVE J0653-56.4, GD 86, EUVE J0823-25.4, 
and EUVE J0916— 19.7. Downward arrows represent upper limits (§3.1). 
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10.00 



0.10 



0.01 

0, 



EUVE J1024-30.3 



10.00 F 



T = 36700 K 



1 1.0 10.0 

X {fj.m) 



0.10 - 



^ 0.01 
100.0 



EUVE J 1439+75.0 



T = 36000 K 



1.0 10.0 
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1.0 10.0 



^ 0.01 L_ 
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Fig. 7.— Spectral energy distribution of EUVE J1024-30.3, EUVE J1439+75.0, EUVE 
J1535— 77.4, and EUVE J1546— 36.7. Downward arrows represent upper limits (§3.1). 
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j, I 
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Fig. 8.— Spectral energy distribution of PG 1609+631, PG 1642+413, PG 1658+440, and 
EUVE J1746— 70.6. Downward arrows represent upper limits (§3.1). 
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Fig. 9. — The sum of the number of point-hke source detections as a function of magnitude 
at each wavelength in the IRAC imaged fields of 16 metal-rich white dwarfs with 600 s 
exposure times (Paper I). 
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Fig. 10. — IRAC 3.6 and 4.5 /im color-magnitude diagrams for all point-like sources in the 
fields of the 6 Hyades white dwarfs observed at both short wavele ngths. Also plotted (stars) 
are representative points in the IRAC T dwarf sequence (Tl— T8; IPatten et al.ll2006l ) at the 
46 pc distance to the open cluster. 



